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Abstract 

Using a 13.7 fb~^ sample of continuum two-jet e"'"e~ qq events collected 
with the CLEO detector, we have searched for correlations between A and A 
particles, specifically in cases where the opening angle between the two par- 
ticles is large and each has momentum >1 GeV/c. Such correlations may 
indicate the presence of baryon number conservation at the primary quark 
level. A previous CLEO study of Ac Ac correlations |^ indicated direct, asso- 
ciated production of primary charmed baryons Ac! e~^e~ cc ^ AcAc. That 
effect was not observed in Monte Carlo simulations. Our current search for 
similar direct, associated production of A baryons at the primary quark level 
{e~^e~ — > ss ^ AA, e.g.) qualitatively indicates a similar effect, although it 
relies on a Monte Carlo dependent subtraction of background AA production 
(based on the default JETSET 7.4 event generator). 
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I. INTRODUCTION 

At ^/s ~10 GeV, below the threshold for e^e~ — > BB, particle production in e~^e~ annihilation 
occurs in a largely low-Q^, non-perturbative regime. Fragmentation models are therefore used to 
describe the process e+e" — >hadrons. Computer codes such as JETSET have been extremely 
successful at matching experimental results on inclusive particle production, multiplicities and 
angular distributions, both qualitatively and, to a large degree, quantitatively. Compensation of 
baryon number is one of the more subtle aspects of particle fragmentation modeling. One obvious 
question is whether baryon compensation occurs locally (e.g., small opening angle between baryon 
and antibaryon) or globally (large opening angles). In the case when a baryon is produced in the first 
step of fragmentation (e"*"e~ — > cc; c — > Ac, e.g.), it is possible that both baryon and flavor quantum 
numbers will be compensated in the opposite hemisphere^ (e.g., e+e" cc; c —>■ Ac; c — > Ac). 
However, in models in which the primary quark and anti-quark fragment entirely independently of 
each other, no such correlations are expected. 

Previous LEP studies (2|-Q of AA production at y/s=90 GeV found that in events containing 
both a A and a A, the two particles tend to be produced with small opening angles between them. 
Those analyses also discrimininated between different models of AA production. It was found in 
the DELPHI and ALEPH analyses that the JETSET string fragmentation event generator gave 
excellent agreement with data when the "popcorn" control parameter p was set to 0.5; the OPAL 
analysis favored slightly higher values of p. (In the default CLEO Monte Carlo event generator, 
we have used p=0.5.) In those previous studies, no statistically significant signal for correlated 
opposite-hemisphere, primary (AA) production was found. 

A previous CLEO study ^ of charmed baryons sought to discriminate between independent 
vs. correlated fragmentation models. Assuming that primary particles fragment independently, 
then the number of times that we find a Ac baryon opposite a Ac antibaryon in an event (denoted 
"Ac I Ac"), scaled to the total number of observed Ac (denoted "(A^AeI")^ should be equal to the 
number of times that we find a Ac baryon opposite any other anti-charmed hadron He scaled 
to the total number of observed anti-charmed hadrons (i^^^). it was found that, given a Ac 
(reconstructed in five different decay modes), a Ac is observed in the opposite hemisphere (0.72 ± 
0.11)% of the time (not corrected for efficiency). By contrast, given a D, a Ac is observed in the 
opposite hemisphere only (0.21 it 0.02)% of the time. Normalized to the total number of either Ac 
or D "tags", that study concluded that it is 3.52ib0.45±0.42 times more likely to find a Ac opposite 



^We define two particles to be "opposite hemisphere" if their opening angle exceeds 90 degrees. 
This definition is therefore disconnected from momentum flow in the remainder of the event and is 
only indirectly related to such standard parameters as thrust, sphericity, etc. As is demonstrated 
in the text, most A baryons emerge either very close to, or directly opposite, the anti-A studied, so 
the separation of an event into hemispheres is at least approximately valid. 
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a Ac than opposite a D meson. This enhancement is not produced in the default JETSET 7.4 
e^e^ cc Monte Carlo simulation. 

As a straightforward extension of that analysis, we can search for similar correlations between 
primary A and A baryons. In this case, the correlation is obscured by the fact that, unlike Ac 
baryons, A's produced in e~^e~ annihilations do not necessarily contain the primary quarks, and 
will be produced copiously in fragmentation, as well as in weak decays of charmed baryons. In our 
current analysis, the production of (A|A) through fragmentation is modeled using the JETSET 7.4 
event generator combined with a GEANT-based Q simulation of our detector. The "feeddown" 
contribution from correlated primary (Ac|Ac) production (Ac — > A|Ac — > A) is evaluated from the 
data itself. 

For this study, A+'s are fully reconstructed in the decay modes A+ — > pK~'ir~^, A+ — > pKg, 
A+ — > Att"*", A+ Avr+vr^vr^", A+ — > pJC^7r"'"7r~,Q and partially reconstructed through A^ AX. 
A baryons are reconstructed in A ^ pTT~ . 

II. APPARATUS AND EVENT SELECTION 

This analysis was performed using the CLEO II and the upgraded CLEO II. V detectors oper- 
ating at the Cornell Electron Storage Ring (CESR) at center-of-mass energies ^/s = 10.52-10.58 
GeV. The event sample used for this measurement is comprised of 9.2 fb~^ of data collected at 
the T(4S) resonance and 4.6 fb"^ of data collected about 60 MeV below the T(4S) resonance. 
Approximately 20 x 10^ continuum cc event^ are included in this sample. 

For 4.6 fb~^ of the data used for this analysis ("CLEO-II" data 0), measurements of charged 
particle momenta were made with three nested coaxial drift chambers consisting of 6, 10, and 51 
layers, respectively. In a subsequent upgrade ( "CLEO-II. V" corresponding to the remaining 
data used for this analysis), the innermost tracking chamber was replaced with a high-precision sili- 
con detector, and the gas in the main tracking volume was changed to provide better cell resolution 
and improved specific ionization (dE/dx) resolution Q. The entire tracking system fills the volume 
from r=3 cm to r=l m, with r the radial coordinate relative to the beam (z) axis. This system is 
very efficient (e >98%) for detecting tracks that have transverse momenta (pt) relative to the beam 
axis greater than 200 MeV/c, and that are contained within the good fiducial volume of the drift 
chamber (| cos 9z\ <0.94, with 9z defined as the polar angle relative to the beam axis). This system 
achieves a momentum resolution of {Sp/p)'^ = (O.OOlSp)^ -|- (0.005)^ {p is the momentum, measured 
in GeV/c). Pulse-height measurements in the main drift chamber provide specific ionization reso- 
lution of 5.0% (CLEO II. V) or 5.5% (CLEO II) for Bhabha events, giving excellent K/tt separation 
for tracks with momenta up to 700 MeV/c and separation of order 2a in the relativistic rise region 
above 2.5 GeV/c. Outside the central tracking chambers are plastic scintillation counters, which 
are used as a fast element in the trigger system and also provide particle identification information 
from time-of- flight measurements. 

Beyond the time-of- flight system is the electromagnetic calorimeter, consisting of 7800 thallium- 
doped Csl crystals. The central "barrel" region of the calorimeter covers about 75% of the solid 
angle and has an energy resolution which is empirically found to follow: 



^Charge conjugation is implicit. 

^Corresponding to approximately 20 x 10^ uu, 5 x 10^ dd, and 5 x 10^ ss events, respectively. 
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(1) 



E is the shower energy in GeV. This parameterization includes noise effects, and translates to an 
energy resolution of about 4% at 100 MeV and 1.2% at 5 GeV. Two end-cap regions of the crystal 
calorimeter extend solid angle coverage to about 95% of 47r, although energy resolution is not as 
good as that of the barrel region. The tracking system, time-of-flight counters, and calorimeter are 
all contained within a superconducting coil operated at 1.5 Tesla. An iron flux return interspersed 
with proportional tubes used for muon detection are located immediately outside the coil and in 
the two end-cap regions. 

Primary proton, kaon or pion charged track candidates must pass the following restrictions: 

(a) The track must pass a 99% probability criterion for its assumed particle identification, based 
on the associated charged track's specific ionization measured in the drift chamber. 

(b) The track must have momentum greater than 100 MeV/c. 

Each reconstructed charmed hadron must have momentum greater than 2.3 GeV/c to ensure 
that there is no contamination from B-meson decays to charm. 

III. DATA ANALYSIS 
A. Production Ratios 

We define the single-tag yield to be the number of reconstructed events containing one particular 
hadron H. This yield is typically determined by fitting a double-Gaussian signal atop a smooth, 
low-order polynomial background function. The number of double-tags is defined as the number 
of events in which two specific particles are both reconstructed, separated by less than 90 degrees 
{HiH2\'-, 'same-hemisphere') or greater than 90 degrees {Hi\H2', 'opposite hemisphere ').Q 

To suppress possible contributions from e~^e~ — > T(4S') — > BB, with subsequent decays such 
as S — > Hc(— > AX)AX, we have imposed a minimum momentum requirement p\ >1.0 GeV/c. 
Decays of S-mesons should generally produce lower momentum A and A; BB events are also likely 
to have different AA production dynamics compared to AA production resulting from direct e~^e~ 
annihilations. (The requirement p\ >1.0 GeV/c is "standard" in other continuum Lambda studies 
as well [p|, [To| , [rT|| .) In order to check whether there is BB contamination of the AA data sample, 
we compare our AA yield derived from data taken on the T(4S) resonance with the AA yield 
obtained using data taken on the continuum below the resonance. We expect these yields to be 
in the ratio of the integrated luminosities of the two samples, corrected for the 1 /s dependence of 
the continuum cross-section, if there is no BB contamination (these two effects give an expected 
ratio of 2.01 for our data). Requiring p\ >1.0 GeV/c, the yield of same-hemisphere AA| pairs on 
the T(4S) resonance (5610.4 it 79.1) compared to the yield for the continuum (2815.2 it 56.1) is, 
indeed, consistent (at the 2a level) with this ratio. Opposite hemisphere yields (4962.8 ± 78.8 and 
2398.9 lb 54.4, respectively) are similarly consistent with an exclusively continuum origin of our A|A 
sample. For maximal statistics, we use our entire sample for the subsequent analysis, and discuss 
residual BB contamination effects later in this document. 



^We use the notation ^^HiH2^ (without a vertical bar) to indicate a generic correlation event in 
which the two particles may be found in either hemisphere relative to each other. 
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FIG. 1. Double-tag invariant mass of Ac candidates plotted vs. invariant mass of Ac candidates 
(Ac|Ac) from data. Shown is the sum of the modes: A+ — > pK^-n^ , A+ — > pK^, A+ Avr^, 
A+ A7r+7r^7r+, and A+ — > pK^n^-K^ (and their charge conjugates, in the case of Ac reconstruc- 
tion). 

Table | summarizes yields in both data and Monte Carlo simulations of comparable size. No 
detection efficiency corrections have been applied. The number of observed inclusive, single-tag 
particles is presented, in addition to the number of same hemisphere double-tags ('AA|'), opposite- 
hemisphere double-tags ('A|A'), and the rate of double-tags per single-tag ('(A|A)/A'). Where 
appropriate, differences of data ratios minus Monte Carlo ratios are given to allow direct compar- 
isons (column 4 in the Table). 

Comparing the Monte Carlo vs. data yields in Table I, we note: 
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FIG. 3. Double-tag plot of A|A (opposite hemisphere) for data. 
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TABLE I. Yields in data vs. Monte Carlo simulations. The data are drawn from a sample of 
~ 50 X 10^ hadronic events (~2/3 on-T(4S) events plus ~l/3 taken on the continuum below the 
T(4S)); the MC is drawn from a sample of 60 x 10^ exclusively continuum events generated by 
JETSET 7.4, fully simulated in our CLEO detector and reconstructed using the same algorithms 
as applied to data. Derivation of final results (specifically, the derivation of the A|A (qq) signal 
yield and evaluation of the non-primary A|A (np) and charmed-baryon A|A (cc) backgrounds) are 
discussed in detail in the text. The difference (Data-Monte Carlo) in the fourth column is shown 
only for the yield ratios, which are sample-size independent. 





Data 


Monte Carlo 


Data - MC 


Ac + Ac {p >2.3 GeV/c) 


83955 ± 1852 


92731 ± 810 




(A + A) (0N-4S, p > 1 GeV/c) 


491087 ± 822 






(A -1- A) (Continuum, p> 1 GeV/c) 


236162 ± 575 


973657 ± 1121 




A + A (total, p> I GeV/c) 


727249 ± 1003 


973657 ± 1121 




(AelAe) 


323.6±39.3 


97.2 ± 40.4 




(Ae|A) + (Ac|A) 


1470.6 ± 74.1 


695.8 ± 67.1 




(AcA|) + (AcA|) 


249.4 ± 26.2 


210.7 ±28.9 




(A|A) 


7361.8 ±95.7 


6519 ±91.0 




(AA|) 


8425.6 ±97.0 


11394.9 ± 112.4 




(AA|) 


28.6 ± 15.3 


36.4 ± 12.4 




(A|A) 


239.8 ± 37.2 


401.0 ± 33.9 




(Ae|Ae)/(Ae + Ac) (xlQ-^) 


3.91 ± 0.47 


1.05 ± 0.44 


2.86 ± 0.64 


(AcA|+AeA|)/(A + A) (xlO-4) 


3.43 ±0.36 


2.16 ±0.30 


1.27 ±0.46 


(Ac|A + Ac|A)/(A + A) (xlO-4) 


20.2 ± 1.0 


7.1 ±0.7 


13.1 ± 1.2 


(A|A)/(A + A) (xlO-2) 


1.01 ±0.01 


0.67 ±0.01 


0.34 ± 0.02 


(AA|)/(A + A) (xlO-2) 


1.16 ±0.01 


1.17 ±0.01 


-0.01 ±0.01 


(AA|)/(A + A) (xlO-5) 


3.93 ±2.10 


3.73 ± 1.27 


0.20 ± 2.45 


(A|A)/(A + A) (xlO-5) 


33.0 ±5.1 


41.2 ±4.1 


-8.2 ±6.6 


Maximum (A A) (np) background 


4820 ± 82 


6519 ±91 




Tagged MC (A|A) from (Ac|Ac) evts. 




477 ± 22 




(A|A) (np) estimate, corrected 


4466 ± 84 


6042 ± 93 




Maximum (A|A) (cc) background 


1671±221 






(A A) (qq), max. bkgnds. 


872±288 






(A|Ac) evts. from (DA)\Ac 


200 ± 34 


169 ± 13 (tagged) 




(A A) (cc), corrected for ((DA) Ac) 


1251±223 






(A A) (qq), (A|A) (cc) corrected 


1290±371 






(A|A) (qq), ALL corrections 


1643 ± 372 
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FIG. 4. Double-tag plot of AA| (same hemisphere) for data. 

1. There is an enhancement, in data, of the number of Ac|Ac per event, relative to Monte Carlo 
(top line of Table I). This correlated production was the subject of our previous paper |5|, 
and was interpreted as evidence for correlated production of charmed baryons from primary 
quarks .0 

2. There is an enhancement in the Ac|A yield in data relative to Monte Carlo. This can largely 
be attributed to the aforementioned correlated Ad Ac production, in which the decay Ac ^ A 
results in a A in the hemisphere opposite the Ac. Such events will also result from cases in 
which charm is compensated by a meson rather than a baryon: h.cK\KD. 

3. In four-baryon events (events with either AA| or A|A), both data and Monte Carlo show a 
preference for A|A over AA|. This is consistent with a model in which two AA| pairs, each 
with cosOpj^ — >1, are created such that one AA| pair is directly opposite the second AA| 
pair. The actual source of these events, although interesting, is not the focus of our current 
effort and its discussion will be deferred until further study. For now, we note that to the 
extent that AAAA events are contributing to both AA and AA samples, statistical consistency 
between Monte Carlo and data in the relative ratios of (AA|)/A and (A|A)/A indicates that 
Monte Carlo simulations model such events reasonably well. 

4. The number of same-hemisphere AA| correlations, divided by the total number of single-tag 
A particles, gives ratios, respectively for data and for Monte Carlo: 



We have loosened the Ac selection requirements for this analysis relative to our previous analysis. 
This results in approximately 20%(44%) larger single-tag(double-tag) Ac reconstruction efficiency. 
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DATA : = (1.16 ± 0.01)% 

(A + A) ' 

MONTE CARLO : = (1.17 ± 0.01)% 

(A + A) ' 

which are in excellent agreement. This agreement gives us confidence that the Monte Carlo 
simulation provides an adequate model of the non-primary component, which is expected to 
dominate the small opening angle sample.^. 

5. By contrast, the opposite hemisphere yields: 

DATA : = (1.01 ± 0.01)% 

(A + A) ' 



MONTE CARLO : ^ ' 1. = (0.67 ± 0.01)% 
(A + A) ' 

indicate large opening-angle A|A production in data at a rate 50% greater than the Monte 
Carlo simulation]^ The total number of observed A|A events should arise primarily from 
three sources: 

a) our signal of interest: direct primary production in e^e^ ^ QQ (q = u, d, or s) events by a 

mechanism similar to that which leads to the observed Ad Ac enhancement (designated 
"(A|A)(qq)"); 

b) A|A events due to feeddown from e~^e~ cc ^ AcAc, with: Ac — > (A|Ac A) (desig- 

nated "A|A(cc)") and 

c) non-primary A baryons which do not contain a primary light quark and are not Ac decay 

products (designated "A|A(np)"). 

Generally, we will use "fragmentation", or "non-primary" to denote particles not containing 
primary quarks and particles which are not direct weak decay products of charmed hadrons. In 
contrast to a) and b), fragmentation will lead to events containing both A and A in the same- 
hemisphere, as well as opposite- hemispheres. We now outline our subtraction of these background 
components. 



^We do not include among the non-primary hadrons the feed-down A decay products of Ac, since 
they include part of the primary hadron. 

was found that setting the JETSET parameter p=0.6, as suggested by OPAL, resulted in 
a larger A|A yield, but also produced a substantial deficit in the same-hemisphere Monte Carlo 
AA| yield compared to data, and thus gave a considerably inferior match to the AA fragmentation 
component compared to p = 0.5. 
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B. Subtraction of non-primary AA component 



To estimate the non-primary AA contribution, we rely heavily on JETSET 7.4 Monte Carlo 
simulations. We use the following procedure to subtract AA (np) production in our search for an 
opposite-hemisphere, correlated, primary A|A (qq) signal: 

1. We plot the cos O^j^ distribution in both data and Monte Carlo simulations, with Oj^-^ defined 
as the opening angle between A and A. 

2. We normalize the Monte Carlo AA yield in the forward hemisphere (cos^^;^ >0) to data, and 
subtract the result from data, for the full (cos d/^-]{) angular distribution. As noted previously, 
the match (absolutely normalized to the total number of A's) between data and Monte Carlo 
simulations is satisfactory in the forward hemisphere (cos^^;^ > 0), however there is an 
under-subtraction of events in the back hemisphere. This subtraction therefore results in an 
excess of back-to-back A|A events in data relative to Monte Carlo simulations. 



C. Subtraction of Ac|Ac feeddown component from A|A yield 

Based on the observed number of: a) Ac|Ac, b) A|A, and c) Ac|A in data, we can calculate the 
total number of A|A(cc) correlations, assuming that any A opposite a Ac is a Ac decay product (we 
discuss below two corrections to these assumptions). This is done by setting the probability that 
both charmed baryons in an e~^e~ cc ^ Ac|Ac decay to lambda baryons, relative to just one 
decaying to a lambda baryon, equal to the probability that one decays to a lambda, relative to the 
probability that neither decays to lambdas. Designating Ac — > A as the inclusive A yield from Ac 
decay, this probability condition can be written as: 

iV(Ac ^ A|Ac A) ^ iV(Ac ^ A|Ac) 

iV(Ac ^ A|Ac) iV(Ac|Ac) ^ ^ 

Monte Carlo simulations indicate that, with the kinematic requirements we impose on our 
candidate event sample, (93 ± 2)% of the final state A's produced in Ac — > A are, indeed, opposite 
each other. Figure ^ shows the cosine of the opening angle between the final state A relative to 
the final state A in such feed-down events; as expected, the distribution peaks at cosS^;^ — 1. 
Thus, assuming that all the observed (A|Ac) events result from e~^e~ AcAc events in which 
Ac — > A in the hemisphere opposite the Ac, the number of "feed-down" (Ac — > A|Ac — > A) events 
can be determined from the above equation directly as: A^^(A|Ac)/A^(Ac|Ac)= 1671 ± 221.0 Using 
this value for the A|A(cc) background and subtracting the scaled MC background as described 
previously, we obtain, using the yields from Table 1: 

NiA\A{qq)) = 

N{A\AUata - N{A\A)mc X ^^^^^ - N{A\A){cc) = 
8425 

7362 - 6519 x 1647 = 872 ± 288 events. 

11395 



^Note that one must use only one charge conjugate in the numerator of this ratio (A^(A|Ac) = 
735.3) to properly compare rates, per reconstructed Ac. 
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FIG. 5. Cosine of opening angle between A and A in e'^e AcAc — > AA events, from Monte 
Carlo simulations. The distribution has been normalized to unity. 

Thus, under the most pessimistic of background assumptions, we obtain a 3cr excess in our 
estimate of the primary, correlated A|A yield. 

Some of the observed A|Ac events will arise from other sources, such as events in which a A 
baryon compensates the baryon number of the Ac {DA\KAc, e.g). There will also be contributions 
from (He — > A) I Ac and (Oc A)|Ac. The number of feed-down (Ac — > A|Ac — > A) calculated 
through the above prescription therefore yields an overestimate of the true feed-down contribution 
and therefore will yield a lower limit on correlated direct A|A production when the feed-down 
component is subtracted. 

We correct our estimate of the AA(cc) background (1671 it 221) with guidance from Monte 
Carlo simulations. In Monte Carlo simulations, for which the parentage of a given particle is 
known, we observe 169 ± 13 A|A events resulting from Ac|AD events. We estimate the number of 
(A|Ac) events in our data sample resulting from Ad AD events (and thus incorrectly attributed to 
(Ac A| Ac A)) by normalizing the number found in the Monte Carlo by the ratio of (AAc-|- AgA) 
(all angles) in data relative to Monte Carlo simulations: 

/A iTnA iV(AcA + AAc)dato ...A iTnA 249.4 ±26.2 .aa^M-iA 

= Ar(AcA + AAc)MC ' ^^^^'^^^^^ = 211.2 ±28.2 " = ^ 

Using this value, we can now re-calculate the background we attribute to A|A(cc) as: (A/'(A|Ac)— 
199.5)^/A^(Ac|Ac) = 1251±223 events. This smaller background estimate implies a larger correlated 
signal yield: 

A\A{qq) = 7362 - [6519 x (8425/11395)] - 1251 = 1290 ± 371 events. 

To better estimate the full effect we have to carry out the additional subtraction described in the 
following section. 
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D. Additional correction due to Ac|Ac production in Monte Carlo simulations. 

Our subtraction of the non-primary component from Monte Carlo simulations has not been 
corrected for known A|A(cc) contributions in the simulation itself. The opposite-hemisphere A|A 
yield from Monte Carlo simulations given in Table | is therefore an overestimate of the non-primary 
component that we subtract out. Since the parent type in simulations is known, this correction can 
be made directly. We tabulate 477 it 21.8 AA(cc) events contributing to our A|A sample in Monte 
Carlo simulations. The calculated yield of primary, correlated A|A events is now: 

A\K{qq) = 7362 - [(6519 - 477) x (8425/11395)] - 1251 = 1643 ± 372. 



E. Summary of Subtraction Procedure 

Figure ^ displays the cos^^;;^ distributions, as well as the subtraction procedure. Overlaid 
with the data cos 6j^j^ distribution is the Monte Carlo simulation. Applying the maximum possible 
A|A(cc) background estimate to the data distribution shown in Fig. (i.e., not correcting for DAcA 
background in data or A|A(cc) in Monte Carlo simulations), we obtain the most conservative back- 
hemisphere excess of 872 it 288 events (statistical errors only). Applying all corrections (as in the 
Figure), this excess is approximately doubled. 

F. Search for correlated A|A excess relative to AK 

In our study of Ac|Ac correlations, we compared Ac|Ac to Ac\D. In an analogous way, we shall 
attempt to compare A|A production relative to A\K production. To normalize properly, we have 
compared the fractional production rate: (A|A)/A relative to {A\K)/K, where the denominator 
designates the total number of detected single-tag A (or K). This technique was used to search for 
evidence of correlated Ac|Ac production in our previous publication. While a D tag always contains 
a primary quark, our K sample includes non-primary kaons, including those that are decay products 
of charmed hadrons.^ Although, in principle, the D\K yield can provide some guidance as to the 
c ^ K rate, there are still considerable uncertainties resulting from the exact admixture of non- 
primary kaons relative to c — > -fC. What must be true, however, is that ss — > A\K events, in which 
the A and K both contain the primary quarks will result in only one flavor correlation; non-primary 
contributions can give rise to both A\K and A\K correlations. Even in the absence of correlated 
primary ss — > A\K production, however, we still expect, because of strangeness conservation, A 
production in association with K (and not K). For this analysis, we restrict ourselves to continuum 
data to exclude possible B KX contributions. Yields are summarized in Table 0. 

We have searched for evidence of correlated production of (A|A) by comparing the ratio of 
(A|A)/A to {A\K)/K] in each case, we can compare with "wrong-flavor" combinations ((A|A)/A 



^There are further complications arising from cases such as: e^e'^ AA(1420); A(1420)— > pK^ . 
In such cases, although observation of the A opposite the would be interpreted as a primary 
baryon-meson correlation, the true underlying event is a primary baryon-antibaryon correlation. 
We can safely neglect such instances as long as the Monte Carlo simulation is producing such events 
at approximately the correct rate. 
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FIG. 6. cos ^^-^ distribution for data (•) with total background overlaid (unshaded histogram). 
The Monte Carlo cosOj^j^ distribution is corrected for the expected (Ac — > A|Ac — > A) contribution 
(normalization obtained from data, with the shape taken from Monte Carlo simulations) , applying 
all corrections cited in the text. That total background distribution has been normalized to the data 
in the region cos0^-^ > 0, then subtracted from the data distribution. The bulk of the background 
histogram is due to A|A (up) production. The component of the total background due exclusively 
to A|A (cc) is shown as the shaded histogram. The final, background-subtracted (all corrections 
applied) data excess is shown in the lower panel. 
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TABLE II. Kaon-Lambda correlation yields in data vs. Monte Carlo simulations. The data 
are drawn from the same sample as for the previous Table. Derivation of final results is discussed 
in detail in the text. No fake subtractions have been applied to the quoted charged kaon yields. 





Data 


Monte Carlo 


K~ fContinuum ONLV) 


325427 ± 570 


1130Q87 ± 1064 


K+ (Continuum, ONLY) 


330010 ± 574 


1149240 ± 1072 




2009 5 ± 50 6 


9426 3 ± 106 




738.8 ± 30.5 


2598.4 ± 55.2 


{h\K-) 


979.6 ± 38.7 


3933.6 ± 75.4 


{kK-\) 


344.8 ± 23.2 


899.4 ± 37.6 


{k\K+)/K+ (xlO-3) 


6.1 ±0.1 


8.2 ±0.1 


{KK+\)/K+ (xlO-3) 


2.2 ±0.1 


2.3 ±0.1 


{h\K')/K~ (xlO-3) 


3.0 ±0.1 


3.5 ±0.1 


{kK-\)/K- (xlO-3) 


1.1 ±0.1 


0.8 ±0.1 


Vk+ (xlO-3) 


3.9 ±0.2 


5.9 ±0.1 


Vk- (xlO-3) 


2.0 ±0.1 


2.7 ±0.1 


S 


8.8 ±0.3 


3.8 ±0.1 



and {K\K)/K), as well as same hemisphere combinations, as estimates of the fake-kaon and non- 
primary components. Contributions to KK^ will arise from: i) correlated, primary production, ii) 
weak decays of charmed hadrons in Ad-D events, iii) compensation of non-primary strange quark 
production, and iv) fake kaons. Contributions to KK~ will arise primarily from non-primary and 
fake kaons. In order to isolate correlated primary production, we therefore define, for both Monte 
Carlo simulations as well as data, the subtracted, normalized fractions: 1^7^+ ^ (A|i^+)/K+ — 
{kK+\)/K+ and V^- = {k\K-)/K- - {kK-\)/K-. We also define a signal ratio S, in analogy 
to our previous publication S = Note that, in constructing this ratio, we subtract 

the contributions to the A|A enhancement from Ac — A|Ac — A to ensure that any observed 
enhancement is not due to the previously measured Ad Ac enhancement. We find S{dato) = 8.77 ± 
0.26, compared with S{MC) = 3.79 ± 0.08 (statistical errors only). This is consistent with an 
enhancement of correlated primary baryon-antibaryon production compared to correlated primary 
baryon-meson production. We note that the (unevaluated) systematic uncertainties in this analysis 
are expected to be considerable, as indicated by the discrepancy in the T>j^+ values between data 
and Monte Carlo simulations (suggesting a greater likelihood of strangeness conservation to occur 
by production of mesons vs. baryons in simulations), as well as the discrepancy in the T)j^^ values 
(suggesting different fake rates in the simulation compared to the data). Because of the large 
systematic errors associated with this exercise, the significance of the difference between S (data) 
and S (Monte Carlo) cannot be simply evaluated on the basis of the statistical errors. This exercise 
is to be viewed only as a check of the primary A|A production enhancement discussed in subsections 
III A through III E. 
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G. Comparison to correlated Ac|Ac production 



In principle, one might hope to compare the yield of correlated A|A production to that of 
Ac|Ac. Quantitatively, one could compare correlated primary baryon production for charmed vs. 
non-charmed baryons, via the ratio: N(Ac|Ac)/e^^/cc relative to: N(A|A)/e^/qq. The efficiency 
factors e include both the reconstruction efficiency in each case, as well as the fraction of the 
momentum spectrum accepted, given the minimum momentum requirements in each case {p > 2.3 
GeV/c for the Ac analysis vs. p > 1.0 GeV/c in the A analysis) What is unknown in the case of 
the A is what fraction of uti vs. dd vs. ss evolve into p|p, n|ii, A|A, A|S^, A\K^p, etc. Nevertheless, 
ignoring such unknowns, and using efficiencies eA and from Monte Carlo simulations (eA ~ 0.091 
and eAc ~ 0.022), we expect to observe ~5300 A|A correlations, scaling from the observed Ac|Ac 
correlations. We observe ~20% of the expected value, consistent with a model in which primary 
baryon-antibaryon correlations in ss events are populated equally by A|S, A|A and with very 
little contribution from e~^e~ — >uu and e~^e~ — >dd. Unfortunately, a) photon-finding (S — > A7 
or S"^ — > pTT^) systematics, b) substantially lower signal-to-noise ratios compared to A ^ Pt^^ , c) 
reduced statistics, and d) the necessity to reconstruct a secondary vertex from a single charged 
track plus neutrals (S+ Ptt*^, e.g.) make the S correlation analyses considerably more difficult 
than the analysis described herein. 

IV. SYSTEMATICS 

Since this analysis relies crucially on the ability of the Monte Carlo simulation to model angular 
correlations in the data, it is important that related kinematic parameters be checked. To verify that 

JETSET models two-particle opening angle (cos 9^ ) distributions adequately, we have compared 

Monte Carlo expectations to the data for the opening angle distribution between oppositely signed 
tracks (Figure ^). For six different minimum momentum requirements, the Monte Carlo simulation 
is observed to model the data quite well. 

To examine the dependence of our result on the minimum momentum requirement, we have 
analyzed the A|A excess (data - background) using different minimum momentum restrictions on 
our A sample. We find that the match between data and simulation, for the normalized same- 
hemisphere yield (AA|)/A, remains excellent for p > 0.5 GeV/c (0.0164 it 0.0001 for data vs. 
0.0168 ± 0.0001 for simulation) as weh as for p>l.2 GeV/c (0.0092 ± 0.0001 for data vs. 0.0091 ± 
0.0001 for simulation). With all corrections applied, opposite- hemisphere excesses are still observed 
for PA >0.5 GeV/c (3306 ± 543 events) and pA >1.2 GeV/c (865 ± 283 events). 

Although the minimum momentum requirement (pA >1.0 GeV/c) should be highly efficient at 
removing backgrounds from B AX, there is still some residual contamination of our single-tag A 
and our double-tag AA sample from T(4S) decays. Inspection of Table I indicates that these may 
represent ~3% corrections to both the single-tag and double-tag A|A samples; our yields indicate 
that the contamination to the same-hemisphere AA| yield is considerably less. Therefore, making 
such corrections explicit would have the effect of modifying our ratios and our corrected yields by, 
at most, ~3%. 



One might argue that since, in the scaled variable (3 = p/m, the two requirements are approxi- 
mately equivalent, the acceptances are also comparable. 
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FIG. 7. cos(0)-| , defined as the opening angle distribution for all oppositely-signed charged 

track pairs, for data (+) vs. JETSET 7.4 Monte Carlo simulations (-). The histograms have been 
normalized to equal areas. 

The analysis would obviously be cleaner if the parent of each A could be unambiguously identi- 
fied. Because it is not possible to definitively distinguish A baryons which contain primary quarks 
from non-primary baryons, the extraction of the correlated, primary A|A signal is inherently Monte 
Carlo dependent. Only in the previous case of Ac|Ac production could one conclusively distinguish 
first-rank from higher-rank baryon production (since higher-rank Ac production, at our energies, 
is expected to be zero). In principle, one might hope to separate primary A production from 
non-primary A production through several techniques; each technique is, however, fraught with its 
own particular difficulties. Measurements such as three-fold KDKc correlations or e^e~ —>■ jAAX 
production, in which the (initial state radiation) photon has sufficiently high energy to exclude 
Ac AX production, can help us assess, e.g., the contribution to A production from charmed 
baryon decays, but cannot define the presence or absence of a correlated primary signal. In the 
CLEO-c era {^/s ~4 GeV), the limited phase space should make this measurement considerably 
simpler. 

V. E\E CORRELATIONS 

In principle, (H|H) pairs may be used to further refine our understanding of primary quark- 
antiquark production and may, e.g., be used as a discriminant between different [|l^ models of 
baryon-antibaryon correlations. Examining Avr" invariant mass combinations, we have observed a 
small, but statistically significant signal for H|H correlations (Figure |8|). Unfortunately, the limited 
signal size is insufficient to attempt to measure primary correlated H|H production. Nevertheless, 
it is of interest to compare the yield of normalized to the total number of charged cascades: 

H+H^/(H+ +H^), relative to the corresponding value for lambda baryons: AA/(A + A). Loosening 
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FIG. 8. Same hemisphere and opposite hemisphere H+H candidate invariant mass distributions 
(ps >0.5 GeV/c). The correlation is most pronounced in the same hemisphere. 

our minimum momentum requirement to 0.5 GeV/c, we find 14537 it 135 total (H+ + H^), an 
opposite-hemisphere yield = 13.0 ib 3.9, and a same hemisphere yield = 21.2 it 

5.1. Correspondingly, we obtain + H^) = (1.5 it 0.4) x 10~^ for same-hemisphere 

cascades, vs. (AA)|/(A + A) = (23 ± 0.4) x 10"^ for same-hemisphere lambdas. For opposite 
hemisphere baryon-antibaryon correlations, the corresponding numbers are (0.9 it 0.3) x 10^^ and 
(28 lb 0.6) X 10~^. In contrast to di-lambda production, di-cascade production favors (albeit with 
small statistics) same- rather than opposite-hemisphere production. Normalized to the total number 
of baryons, the di-cascade rate (integrated over all angles) is apparently suppressed relative to the di- 
lambda rate. This is consistent with a model in which cascade and lambda production is dominated 
by light quark popping; in such a picture, E~ suppression is therefore a direct consequence of 
strangeness suppression. 

VI. SUMMARY 

Under conservative assumptions, we observe a ~ 3o" (872 it 288) excess of opposite-hemisphere 
A|A production in data compared to the expectations of the JETSET 7.4 event generator combined 
with the full simulation of our detector, and after accounting for feeddown production of AA from 
charmed baryons. With appropriate corrections applied, this excess increases to (1643ib372) events. 
These results are consistent with enhanced correlated, primary A|A production, of the type observed 
previously in Ac|Ac correlations. However, we stress the inherent Monte Carlo dependence of this 
conclusion (not present in the Ad Ac correlation analysis ||5|), and that the complete parameter 
space of the event generator has not been fully explored. Data-taking planned for CLEO-c at D\D 
threshold should be able to more definitively measure such A|A correlations, in a considerably less 
Monte Carlo-dependent manner. 
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